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Abstract 

Given two positive integers n, r, we define the Gaudin function of 
level r to be quotient of the numerator of 



det(^{{xi - yj){xi - tyj) ■ ■ ■ (xi - fyj)) 



i,j=l...n 



by the two Vandermonde in x and y. We show that it can be char- 
acterized by specializing the Xj into the yj variables, multiplied by 
powers of t. This allows us to obtain the Gaudin function of level 1 
(due to Korepin and Izergin) as the image of a resultant under the 
the Euler-Poincare characteristics of the flag manifold. As a corol- 
lary, we recover a result of Warnaar about the generating function of 
Macdonald polynomials. 

ooo o^ ooo o^ o^ o^ ooo 

oooooooooooooooooooooooooooooooo 



1 Gaudin functions of arbitrary level 

Let X = {xi, . . . , Xn}, y = {yi, ■ ■ ■ , Vn} be two sets of indeterminates of the 
same cardinality n. The Cauchy determinant det((a; — ygy pl^Y^ a 

central role in the theory of symmetric functions [Ti] . 

Generalizations of this determinant appear in the calculation of correla- 
tion functions of different physical models. Gaudin [U Ch.IV] obtained the 
determinant 

det{{x-y)-\x-y + ^)-')^^^^^^^, 

7 a parameter, for a bose gas in one dimension. Izergin and Korepin [6] 
solved the Heisenberg XXZ-antiferromagnetic model with the help of 

det{{x-y)-\x-tyr')^^^^^^^, 
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and Kirillov and Smirnov [TJ Th.l] wrote more general determinants. These 
different determinants have led to an abundant literature, in connection with 
different statistical models and combinatorial enumerations, for example the 
enumeration of alternating sign matrices displaying some symmetries [15]. 

In this text, we shall be interested in a purely algebraic generalization, 
having no physical interpretation to offer. With x and y as above, let t be 
an extra indeterminate and r be a positive integer. 

We propose to study the determinant 

det 



{x-y){x-ty)---{x-ry) 
which can 
A(x)A(y) 



xex,yey 

This is a rational function, which can be written 



/?(x,y(l + --- + r)) 



^n(x,y), (1) 



where -F^'(x, y) is a polynomial symmetrical in both xi, . . . , x„ and yi, . . . ,yn, 
A(x) being the Vandermonde nj<i(-^i ~ ^«)' resultant i?(x, y), for two 
finite sets of indeterminates x, y, being the product of differences 
rixex Uyeyi^ - y), and y(l + ■ ■ ■ + f ) being {fyj : 1 < j < n, < z < r}. 

We shall call the function F^(x, y) the Gaudin function of level r, in 
recognition of the pioneering work of Gaudin so well illustrated in his thesis 
[5] and his book [1]. 

Let us introduce n sets of indeterminates 



■= {Vj =yj,yj,---,yj} , j 



,n. 



Recall that, for any pair (a, b) of indeterminates, the divided difference 
relative to a, b is the operator (denoted on the right) 



f{a, b) f{a, b)da,b = {f{a, b) - f{b, a)) (a - b) 



-1 



The determinant M„ = \R{xi,y^) ^| is the image of the Cauchy determi- 
nant under the product of divided differences 



where is relative to the pair yj, y*^^. 



r— 1 5 
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Indeed, for any j, any x, one has 



/ 0\-l '^0 / 0\-l/ ^1 I 0\-l/ 2\ 

[x-y^] — ^(x-i/j) [x-y^] — >{x-y^) [x-y^] {x-y^} 



2\~l 



so that the divided differences relative to the alphabet y-' transform all the 
terms {xi — y'j)'^, inside the Cauchy determinant, into {xi — y'j)'^ ■ ■ ■ {xi — 

yd-'- 

On the other hand, it is classical, and immediate, that the Cauchy deter- 
minant be equal to A(x)A(y)/_R(x, y). It can also be expressed as 



Afx) 



1 



y 



y 



n — 1 



y=yi,---,yn 



Its product by YliRi'^^y^) can therefore be written as 



Afx) 



iyyR{^,r-y^] 



i=l,...,n,j=0,...,n—l 



(2) 



Since each factor i?(x, y*) commutes with any divided difference dl, k = 
0, . . . , r-1, j = 1, . . . , n, the image of ([2]) under Hi ' ' ' ^r-i equal to 
M„n,^(x,yO: 



Afx) 



n i?(x, r) = {yy i?(x, y' - y,)5o ■ ■ ■ dl^i 



(3) 



We shall now have recourse to symmetric functions, referring to the last 
section for more details, as well as [121 [13] • For any non negative integer j, 
one has 

(yy i?(x, y^ - ■ ■ ■ dl^, = S,,oir; y^ - x) , 
where □ = \ n-l, . . . , n-lj = {n-iy , and Sj-ci is a multi-Schur function. 

r 

Writing y* = (y* — x) + x, and expanding by linearity, then 



k=0 



This identity allows us to transform ([3]) into 



Afx 



(4) 



where now the entries of row i are Schur functions in the difference y* — x. 

Going back to the original variables, that is, specializing each y* into 
yi + tyi + - ■ ■ + t^yi, we obtain the following expression of the Gaudin function 
of level r. 



Theorem 1 The function 

/?(x,y(l+--.+r)) 



is equal to 



A(x)A(y) 



det 



{xi - yj){xi -tyj)--- {xi - f'yj) 



A(y) 



S,,u{y^{l+■■■+f)-^) 



i=l,...,n, j=0,...,n— 1 



The original case of Izergin, Korepin is for level 1, and reads 
1 



^n(x,y) 



Sn-i{yi{l+t) - x). 



A(y) 

S'i,„_i(?/i(l+t) 



0, . . . ,5'„_i,„_i(?/i(l+t) -x) 



i=l....,n 



(5) 



For example, for r = 1, n = 3, the function is 
i^'(x,y)A(y) 



So22{yi+tyi - x) Si22{yi+tyi - x) 
So22{y2+ty2 - x) Si22{y2+ty2 - x) 
So22 (ys+tys - x) 5*122 ivs+ty-s - x) 



S222{yi+tyi -x) 
S222{y2+ty2 -x) 
S222iy3+ty3 -x) 



We have given another expression in [TU], separating the variables x and y. 

The determinant det ((tXj — yjj)~^{t'^Xi — yj-^)~^) specializes into 
det ((xj — yj){xf — yj)^^) when t = exp{2n\^^/3). In that case, F^(x, y) 
becomes the Schur function in the union of x and y of index [0, 0, 1, 1, 2, 2, ... , n-1, n-1] 
(cf. [T71 [H]). More generally, the Gaudin function of level r, when r is 



odd, displays such a global symmetry. In that case, {tx — yt 



yt-^-')-' 

is equal to 



{x—y){x^~^'^—y^~^^) ^, and the determinant det ((xj 
A(x)A(y) 

;rp2 '^0,0,/3,/3,-,(n-l)/3,(n-l)/3lX + yj . 



r+1 



X ■ 



U^:^'-yj 

More general determinants displaying a symmetry in xi, . . . ,Xn,yi, 
are given in Lemma 13] and [15]. 



■,yn 
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We shall now characterize the Gaudin function by specialization. Ex- 
panding the determinant expressing F^(x, y) by linearity in x, one sees that 

F^(x, y) is a linear combination of terms y^^, ....... y^" , 

< fi, . . . ,fn < (?T,-l)(r+l), i.e. is a linear combination of Schur functions 

of y indexed by partitions contained in ffl = [ (n-l)r, . . . , (n-l)r]. By sym- 

^ . ' 

n 

metry x ^ y, F.J^{'x, y) is a linear combination of products of Schur function 
of X and of Schur functions of y indexed by partitions contained in ffl. 

Given any infinite set of indeterminates z, any linear combination of Schur 
functions in x with coefficients in z, indexed by partitions C ffl, is also a lin- 
ear combination of Grassmannian Schubert polynomials y„(x, z), v C ffl. As 
such, it is determined by the specializations x C {zi, . . . , Znr+n-r}- 

In the next theorem, we choose z = {yi, . . . ,yn, tyi, . . . , t'^yi, . . . , t'^yn, ■ ■ ■} 
to get simple specializations. 

Theorem 2 F^(x, y) is the only linear combination of Schur functions in 
X, with coefficients in y, indexed by partitions contained in ffl, which has the 
same specializations 

X C {yi, . . . ,1/n, • • • ,^"2/1, • • ■,t''yn} 

than the function 

G;(x, y) := ^ fj Saiy^iU . ■ ■ +f ) - x) , 

where □ is, as before, equal to {n-iy . 

Proof. If the specialization of x contains several occurrences of the same 
yi (ignoring the powers of t), then all the functions Sj^\j{yi{l+ ■ ■ ■ +t^) — x) 
vanish. Thus, F^(x, y) as well as (j5^(x, y) vanish in that case. 

Let now x' = {x[ = yit"^, ■ ■ ■ ,x'^ = ynt""}, with < ei, . . . , e„ < r. In 
that case, each yi{l+- ■ ■+t'^) — x' is equal to the difference of two sets of 
respective cardinalities r, n — 1, and, according to (fTTj) . 



5',, a{yi{l+ ■ ■ ■ +f) - x') = (x: - x') Su{y^{l+ ■ ■ ■ +f) - x) . 

This factorization allows to extract from the determinant expressing F^lx.', y) 
the factor Yli Sn{yi{^+ ■ ■ ■ +t^) — x'). There remains det[Sj{x[ — x')), which 
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is equal to the Vandermonde A(x'). Therefore, the two functions F^(x, y), 
G^{'K,'y) have the same speciahzations in . . . , fy^}. To characterize 
-F^(x, y), we need only specialize x to a subset of the first nr + n — r elements 
of {yi, . . . , t^i/n}, and this finishes the proof of the theorem. QED 

Notice that all the specializations occurring in Theorem [2] are either 0, or 
products of factors {i/it'' — Xj). 

When r = 1, Theorem [2] claims that -F^(x, y) has the same specializations 
X C {yi, . ..,yn, tyi, tyn] as A(x) A(y)"^ Hi Sn-i{yi + % - x). We have 
moreover remarked that one can suppress one letter from . . . ,t?/n} to 
characterize the Gaudin function. 

For n = 2, r = 2, the expression in Theorem [T] specializes, for xi = 
yi,X2 = t'^y2, into 



1 



5oii (tyi + t^yi - t^y2) Sni (tyi + t^yi - t^y2 
Sou (2/2 + ty2 - yi) 5*111 (2/2 + ty2 - yi) 

1 -t'y 
1 -yi 



2/2 - yi 
= t%y2-yi)iy 



ty2 



i (2/2 -2/1) (2/1-^2/2) (r2/2- 2/1) • 



Theorem [2] gives on the other hand 
i^2/2 - 2/1 



2/2 - 2/1 



Souityi + t 2/1 - ^ 2/2) 5'oii(2/2 + ty2 



yi] 



t'2/2 



2/1 



2/2 - 2/1 



{tyi - t^y2){t^yi - ^^2/2) (2/2 - 2/0(^2/2 - 2/1) 



which is, indeed, equal. 



2 Euler-Poincare characteristics 

We go back to the original Gaudin- Izergin-Korepin determinant [U [6] , that 
is, from now on, we take level r = 1. 

The Euler-Poincare characteristics for a flag manifold under G/(C"'), con- 
veniently generalized by Hirzebruch [HI [9], can be combinatorially interpreted 
as a summation over the symmetric group : 

/ rr {tx -x')\'" 
C[t][xi,...,x„]9/^/luj_= ^ Get)m(x). 
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This morphism is characterized by the fact that the images of dominant 
monomials 

:= x^' ■ ■ -x^" , Ai > • • ■ > A„ > 

are the Hall-Littlewood polynomials [H] ca-Pa(x, t), the normalization con- 
stants c\, writing A = 0™''!'"^ ■ ■ ■n™", being 

n nii 

The elementary operators (case n = 2) are 

:= lUis, : / — > f^i = fitXi - Xi+i) di 

and generate the Hecke algebra of the symmetric group (as an algebra of 
operators on polynomials. This is the description of the affine Hecke algebra 
that I used with M.P. Schiitzenberger in [5]). The usual generators Tj : = 
— 1 satisfy the braid relations, and the Hecke relation (Tj — t){Ti + 1) = 0, 
while = (1 + t)iyJi. 

We shall also need an affine operation 6, which is the incrementation of 
indices on the x-variables : 

Xi9 = Xi+i , periodicity a;j+„ = Xit~^ . 

Notice that to define Macdonald's polynomials (see next section), one uses 
the periodicity = qxi with q independent of t. 

Theorem 3 Let x, y he two alphabets of cardinality n, f be a function of a 
single variable. Then 

fix,) i?(x - xi , y)(l - te) ■ ■ ■ (1 - r-^e)iui. 

= ifix,)x2 ■ ■ ■ Xnd, ■ ■ ■ ) F„i(x, y) [n]\ . (6) 

Proof. The LHS, as a function of y, belongs to the space generated by the 
Schur functions of index contained in the partition (ra-l)". It therefore can 
be determined by computing all the specializations 

y C {Xi, . . . , Xn, X2t ^,...Xnt ^} 

(we do not need to take Xit^^). 
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To lighten notations, let us take n = 4. The function 

fix,) R{x2+xs+x,, y)(i - te){i - t^e){i - t'e) 

is equal to 



f{xi) R{x2+X3+X4,y) -t 



f{x2) R{xs+X4+xi/t,y) 

/(Xg) R{x4+Xi/t+X2/t, y) - t^f{x4) R{Xi / t+X2 / t+x^ / 1 , y) 

The sum under the symmetric group can be written 

2^ (^/(Xi)i?(x2+X3+X4,y)^^Y^^ J 



■3 



^ f -1 ^A^(2134)^"' 

2^ (^/(xi)i?(x3+a;4+X2t \y) ^^2134) J 

We shall identify the coefficients of /(xi) in both members of O Consider 
all the specializations y C {xi, . . . , X4, X2/t, . . . , a;4/t} of the LHS. Up to 
symmetry, the only non-zero specializations are 



Af(3214) ^ 
A(3214) 



y - {xi,X2,X3,X4}: -t^i?(2^^±f±^, y) E.:.,=i 

y - {Xi, X2, X3, f }: t^' g] + X4, y) E.:.,=l,.,=4 ( 

y - {Xi,X2, f , f }: -t + X3+X4,y) E.:.,=l,.,=2 (^^^ 

y - {Xl, f , f , f }: i?(x2+X3+X4, y) E.:.,=l 



Up to the global factor 

X2X3X4(l-t)(l-t^)(l-t^) 



At(1234) At(4321)i?(xi, X2+X3+X4) 
8 



these specializations are respectively equal to 

R{Xi+X2+X3,X4) R{Xi+X2, X3+X4) R{xi, X2+X3+X4) 



R{xi+X2+Xs^tX4)' R{xi+X2,tXs+tX4)' R{xi,tX2+tXs+tX4) 

They coincide with the specializations of the RHS of ([6]), thanks to Theorem 
El writing f {xi)x2XsX4did2d3 as 

f{xi)x2X3X4 ^ f{x2)xiX3X4 _^ f{x3)xiX2X4 ^ f{x4)xiX2X3 



R{xi, X2+X3+X4) R{X2, X1+X3+X4) R{X3, X1+X2+X4) R{X4, X1+X2+X3) 

In final, we have checked enough specializations to prove ([6]). Q.E.D. 

3 Generating functions of Macdonald poly- 
nomials 

The symmetric Macdonald polynomials Pa(x; q, t) satisfy a Cauchy formula : 

^ ^ ^ xex,j/Gy i>0 ^ ^ A 

sum over all partitions of length ^{X) < n, the constants b\ being defined in 
[a VI.4.11]. 

Let Tq be the following incrementation of indices on the x-variables : 

Xi Tq = Xi+i , periodicity Xi+„ = qXi . 
We want to compute 

Since 

1-t 1-t 

X- Tq = X- h Xi(t - 1) , 

I — q ^ — q 

one has 

ai ( ) (1 - tTq) ■■■{!- rTq) 



l-q 



ai (xy(l-t)) (1 - tro) ■ ■ ■ (1 - tVo) a, [ xyg^^^) (8) 
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The parameter q has been ehminated from the operation, and we are thus 
reduced to the case of Hall-Littlewood polynomials, which is treated in the 
next theorem. 

Theorem 4 The image of the generating function of Hall-Littlewood poly- 
nomials Pxi^, t) = Pa(x; 0, t) under (1 — tro) ■ ■ ■ (1 — t"ro)iyJ^ is 

ai (xy(l-t)) (1 - tro) ■■■{!- t^To)U^ = ai(xy) F„i(x,y) [n]\ , (9) 

where i^^(x, y) is the Gaudin function (xi ■ ■ ■ F^(x^, y), 

= {x^\...,x-^}, and [n]\ = (1 - t) ■ ■ ■ (1 - T). 

Proof. One rewrites (Xi (xy(l-t)) = /2(tx, y^)i?(x, y^)^^. Notice that 

RjtXj^ htxn,y^) _ RjtXj H h tXn, y^) RjtXj+i H h y^) 

= {R{tXi,y'')-t R{xi,y ))—^. ; ■ ^ = ^7 ; ■ ^ 

^ R{xi^ hx„,yV) R{xi^ hx„,yV) 

with /(xj) a polynomial in Xj of degree n — 1. Multiplying the LHS of ([9]) by 
the function i?(x, y^), one transforms it into 

{j{xi)R{tx2 H h tx„, y^) - eif{x2)R{txz H h te„ + xi, y"^) 

+ ■ ■ ■ + (-l)""ie„_i/(x„)i?(xi + ■ ■ ■ + x„„i, y^; 

where Ci, . . . , e„_i are the elementary symmetric functions oft, ... , t^~^. One 
recognizes in this last expression 

f{xi)R{tx2 + ■■■ + tXn, y^)(l -t9) ■■■{!- f'-'e) . 

We now invoke Theorem [61 Since the function /(xi)x2 ■ • • Xndi ■ ■ ■ dn-i is 
a constant (for degree reasons), the LHS of (IHl) is proportional to F^(x, y^), 
that is, is proportional to F^(x, y). In fact, /(x) = Hjl^^ ~ Vj^) 
-t"Y{j{x-y-^), so that 

/(Xi)x2 ■ --Xudi ■ --dn-l = (1 - r)(?/i • --yny^ ■ 

Correcting by the right powers of (xi ■ ■ -Xn) and {yi ■ ■ -yn), one finishes 
the proof of the theorem. Q.E.D. 

One can now go back to the case of Macdonald polynomials, and recover 
a result of Warnaar [TSl Th.3.1]. 
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Theorem 5 (Warnaar) There holds 

A i=l \ 1/ 

(10) 

Proof. The non symmetric Macdonald polynomials are eigenf unctions of cer- 
tain commuting Dunkl-type operators ^i, . . . ,^n, first introduced in [T] and 
extensively used by Cherednik [2]. The eigenvalues are q^^t"'~^, . . . , q^"t^ for 
the polynomial Mx indexed by A : Ai > ■ ■ ■ > A„ > 0. Up to normalization, 
the image of Mx under lyj^ is equal to Pa(x, g, t), and for any symmetric 
function g in n variables, then 

Pa(x; g, t) g{^u ...,U = ^a(x; g, t) g{q^H''-\ . . . , q^H^) . 

Using [11] that lyj^ = lyj^ t"-Vg y^, one sees that 



Pa(x; g, t) - q'-t^-'"-') = ^a(x; g, t){l - tr^) • ■ ■ (1 - tV, 



i=l 



Therefore, the LHS of flTUl) can be identified with 

^1 f^y^l (l-tr,)---(l-rr,) 



1-g 

Thanks to (18]) and (19]), this can be written 



ai(xy)(7i ^xyg^^ ^n(x,y) = (Ti (^y^-!^ ai(txy)P„^(x, y) , 
which is Warnaar's formula. Q.E.D. 



4 Note: Symmetric functions and Schubert 
polynomials 

We use A-ring conventions to describe symmetric functions. Given three sets 
A, B, C of indeterminates ("alphabets"), the generating function of complete 
functions Sn{AB — C) is 

a. := ^ = V z-SMB - C) . 
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We write alphabets as sums of the letters composing them. For example, 

A(l + t H h to is the alphabet {af : a e A,0 <i <r}. 

Schur functions Sv{A — C), v G N", are determinants of complete func- 
tions: 

SM -C) = det(5.^+,_,(A - „ . 

One generalizes Schur functions to multi-Schur functions by taking dif- 
ferent alphabets in blocks of columns of the preceding determinant. For 
example, 5'^i;„2, (^i ~ Ci, A — C) is the determinant with first column 
S^^{Ai — Ci), . . . , S'^j_„+i(y4i — Ci), and entries Sy.+j^i{A — C) elsewhere. 

Multi-Schur functions satisfy some factorization properties [T^ Prop. 
1.4.3]. We need only the following case, which was much used in classical 
elimination theory in the 19th century. 

Lemma 6 Given two finite alphabets A,B of respective cardinalities a,l3, 
given j : < j < /3, then 

S,MA-B) = S,{-B) n ia-b) = {-iye,{B) J] (H) 

aeA,beB aeA,beB 

where ej{B) is the elementary symmetric function of degree j in B. In par- 
ticular, Sj^fic,{A -B) = zfAnB^^. 

There are several families of non-symmetric polynomials extending the 
basis of Schur functions. Of special interest are the Schubert polynomials 
yt,(x, y), V G N", which constitute a linear basis of the ring of polynomials 
in X = {xi, . . . , Xn}, with coefficients in ?/2, . . . , ?/oo- They can be charac- 
terized by vanishing properties related to the Bruhat order. 

The subfamily of Schubert polynomials indexed by (increasing) parti- 
tions form a basis of the ring of symmetric polynomials [12]. It satisfies 
the following property. Given u G N" : < Ui < ■ ■ ■ < Un, let y^"^ : = 

Lemma 7 Let v = [0 < vi < ■ ■ ■ < Vn] be a partition. Then yy(x, y) is 
the only symmetric function in x of degree \v\ = f i + ■ ■ ■ + f « such that 
^^(y<"^,y) = for all u: \u\ <\v\,u=^v, and Y^{x, {0,0, .. .}) = 5^(x). 

More generally, Y^^y^"^-^ ,-y) ^ iff the diagram ofu contains the diagram 
ofv. 

In the case where y = {0, 1,2,.. .} (resp. y = {g", q^, q^, ■■■}), the poly- 
nomials yu(x, y) are called factorial Schur functions (resp. q-factorial Schur 
functions, and the above vanishing properties are extensively used in [TB] . 
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